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Abstract 

We study the gravitational collapse of the axion-dilaton system suggested by 
type IIB string theory in dimensions ranging from four to ten. We extend 
previous analysis concerning the role played by the global SL{2, R) symmetry 
and we evaluate the Choptuik exponents in the elliptic case. 



^ On leave from Department of Physics, Ferdowsi University of Mashhad, P.O. Box 1436, Mashhad, 
Iran. 



1 Introduction 



It is nearly twenty years since M. Choptuik discovered critical phenomena in gravitational 
collapse, and revolutionised the field of numerical relativity [HE]- In this letter we revisit 
the self-similar, spherical gravitational collapse studied in the past in four- dimensions in 
[21 H] , and extend their analysis to higher dimensions and to other implementations of self- 
similarity. Our work is motivated by the AdS/CFT correspondence [3 El (TJ IH] ■ We are still 
far from a holographic description of black hole formation in type IIB string theory. Some 
exploratory work was done in (Qj [10] . In the context of type IIB (for reviews of string theory 
and duality see [HI [12]) and AdS/CFT one would like to consider collapse on spaces which 
approach asymptotically AdS^ x S^. With the bosonic fields in the theory, a natural system 
to consider involves the axion-dilaton and the self-dual 5-form field. One can show that in 
five dimensions the simplest dynamical set to study involves just the Einstein- axion-dilaton 
system with a cosmological constant. This poses an apparent problem in considering self- 
similar collapse, because Einstein spaces do not admit homothetic vector fields. However, 
in the context of critical gravitational collapse we are considering the collapse of matter to 
form small mass black holes and we only need to consider a small space-time region close to 
where the singularity forms. This should be independent of the asymptotic structure of the 
space-time where the collapse takes place. There is numerical evidence in asymptotically 
AdS space-times showing that this is the case [13] . Hence we will eliminate the cosmological 
constant and analyse self-similar critical collapse in dimensions 4 — 10. We are interested in 
the critical solution and the corresponding Choptuik exponent. Our analysis follows closely 
[31 H], and the physical and geometrical interpretation is essentially the same, so we will 
not review their work in this short letter. 

In section two we write the action and the equations of motion of the axion-dilaton 
system, and the conditions that follow from the requirement of continuous self-similarity. 
We found three possible conditions and associated ansatze for the matter fields. We called 
them the elliptic, hyperbolic and parabolic cases associated to the three different clases of 
SL{2, R) transformations that can be used to compensate for a scaling transformation in 
space-time. The analysis in [3l H] corresponds to the four-dimensional elliptic case, but 
their methods can easily be adapted to the two other CcLSGS cLS well as to other dimensions. 
We do not strive to obtain very precise numerical results given that we have no practical 
applications for them. We are interested in exploring the existence of critical solutions and 
the values of their Choptuik exponents as a function of the space-time dimension. For a 
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massless scalar field this has been done in [HI [15], the shape of the plot of our exponents as 
a function of dimension is quite analogous to the results of Oren and Sorkin in [15], however 
more precise numerical work would be necessary before a detailed comparison is possible. 
More details as well as the results for the hyperbolic and parabolic cases will appear in a 
subsequent publication [16j . 



2 The axion/dilaton system 



We will consider the d-dimensional spherical collapse of an axion-dilaton (a, $) system. 
Both field can be combined into a single complex scalar field r = a + ie~^ . The action for 
the model is: 

IGttG J ^ ^ y 2(Imr)2y ^ ^ 

where R is the scalar curvature. The equations of motion are: 

Rab- A 9 i^ardbf + dafdbT) = (2) 

V»V„r + !^I^ = 0, (3) 

Imr 

The theory is classically invariant under SL{2, R) transformations: 

ar + b 

where (a, b, c, d) E R, ad — be = 1 and gab does not transform. This group is supposed to 
be broken to its integer subgroup SL{2, Z) as a consequence of electric- magnetic duality in 
String Theory [IIIIIS]. 

Following [3] the spherically symmetric metric can be written as: 

ds^ = (1 + u{t, r)) {-b{t, r fde + dr^) + r^dnl_^ . (5) 

The time coordinate is chosen so that spherical collapse on the time axis first occurs at 
t = 0, hence the collapse takes place for t < 0. We can still implement time redefinitions 
for (E]), hence as in [3j we set 0) = l,t < and regularity for t < implies that 
M(t,0) = 0,t < 0. 

Continuous self-similarity means the existence of a homothetic Killing vector ^ gener- 
ating global scale transformations: 

C^^Qab = "igab ■ (6) 
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In spherical coordinates ^ = td/dt + r d/dr. Defining the scale invariant variable z = —r/t 
self-similarity of the metric means that the unknown functions u{t , r) , b{t , r) are just func- 
tions of z. The next question to address is the transformation of r(t, r) under scale transfor- 
mations. Since the action is SL{2, i?)-invariant, we can compensate a scale transformation 
of the coordinates {t,r) by an SL{2,R) transformation. Thus if we change variables to 
{t, z), we obtain a differential condition for T{t, z): 

d 

t — T{t,z) = CKo + CKi T Q;2 (7) 

with q;o,i,2 real numbers. The quadratic polynomial on the right hand side has two roots 
which can be two complex conjugate numbers, two distinct real numbers or a double real 
root. They correspond to compensating the scaling transformation with respectively an 
elliptic, a hyperbolic, or a parabolic transformation in SL{2,R). By straightforward ma- 
nipulations and redefinitions the three cases yield the following ansatze for r(t, z). In the 
eUiptic case: 

, , .1- (-ty^fiz) 

T(t, r) = I- — . , ( , 8 

under a scaling transformation t ^ Xt, T{t, r) changes by a SL{2, R) rotation. 
In the hyperbolic case: 

r(t r) = ' - ^-'^^^^'^ (9) 

and under a scaling transformation i — >■ At, T{t,r) changes by a SL{2,R) boost. Using a 
SL{2, i?)-transformation we can transform T(t, r) — (— t)'^ fi^)- Finally, in the parabolic 
case a scaling transformation can be compensated by a translation. After some simple 
SL{2,R) transformations of T(t,r) the parabolic ansatz becomes: 

T{t,r) ^ f{z)+ujlog{-t) ; weR (10) 

We have written —t throughout the previous equations because the collapse will take 
place for t < 0, we can also extend the solutions for i > 0, hence the correct ansatze for 
any value of t is to replace —t^\t\. In all these expressions, f{z) is an arbitrary complex 
function and ou a arbitrary real parameter to be fixed by requiring the critical solution to 
be regular. 

3 Equations of Motion 

We collect in this section the equations of motion for the three ansatze, although later we 
will only present results for the elliptic case. Since we are implementing spherical symmetry. 
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the gravitational degrees of freedom do not propagate, there are no gravitational waves. 
Hence u{z),b{z) should be expressed in terms of f{z). Common to all three cases, it is 
simple to show that 

= - {d-3)bizy ^''^ 

This follows from the Einstein equation for the angular variables. The other equations 
of motion for the self-similar solution involve b{z),f{z). In fact the equation for 6 is a 
first order linear inhomogeneous equation for 6(2;)^ whose initial condition is determined by 
6(0) = 1 together with the initial conditions for f{z), f'{z), that are determined by requiring 
smoothness of the critical solution. This determines not only the initial conditions but also 
the possible value of w. 

The equations of motion are quite complicated, and their derivation is rather tedious 
but straightforward. The equations in the elliptic case are given by: 

= ^' + - rr-^m^Mf - If) - ^^^^'^'^ 



(^_2)6(-l + |/|2)2^ ^ (rf_2)6(-l + |/|2)2v^^ ' (rf_2)6(-l + |/P)2' 

_ 2z{b' + z') 2 p , 2 / _ iu{b^ + z^) \ - 2 

^ (^_2)62(_l+|/|2)2^ ^^(1_|/|2)1^ (^_2)62(1_|/|2);^^ 



2iuj{b^ + 2z^) ^ ^, 2 f{z^- ^^^) icoz\l + I/I 



2^ 



^(ci-2)62(-l+ |/|2)2-^-^/ + ^ (Z2_52) + (52_^2)(1_ |J|2) 

I 2a;V|/P \ 2ujh 



(rf - 2)62(62 -22)(1_ I J|2)2y ^ (rf-2)62(-l + |/|2)2 
2iuj [I ia;(l + |/|2) a;2z2|/|2 



(62 - ^2) \^2 2(1 - |/|2) {d - 2)62(-l + |/|2)2 



/• (12) 



It is useful to note that the equations are invariant under a global redefinition of the phase 
of f{z). We can choose the phase of f{z) to any convenient value at a particular z that we 
will choose to be the origin. 

In the hyperbolic case, the equations of motion are quite similar, but now they are 
invariant under a constant scaling / A/, hence we can choose |/(-2)| or its real or 



4 



imaginary part as we wish at a particular value of z\ 



(d - 2)b(f -ff-' ' (d- 2)b{f - SY " ' (d- 2)6(/ - /) 



2 



{d-2)h\f-fr if-f)\f {d-2)h\f-f) 



(rf_2)62(/_/)2- ^ (^2_^2) (52_^2)(J_J) 



{d - 2)62(62 _ ^2)(j _ j)2y ^ _ 2)52(/ _ /)2 



/• (13) 



(62-^2) ^ 2 2(/-/) (rf-2)62(/_/)2 
Finally in the parabolic ansatz the equations of motion are: 



[d-2)hU-f?'' (rf-2)6(/-/)2- (d-2)6(/-/)2 
n = 2z{h^ + z^) 2 ( ujh^ + z^) \ 

^ {d-2)h^U-ff^ ^^{f-f)\ ^ ^d-2)h\f-f))^ 
2u{h^ + 2z-) 2 ( (z^-^) 2C..2 



(ci-2)62(/_/)2^ ^ ,y (^2_52) ' (52_^2)(J_J) 

/ _7 o\ 79/79 ^N/i- i'\9 ) / 7 r>\7_9//' f\9>' 



{d - 2)62(62 _ ^2)(j _ j)2y ^ _ 2)62(/ - fy 

2ijJ ( \ CO LiP'Z^ 



(b^-z^) V 2 (/-/) (rf-2)62(/-/)^ 



(14) 



They are invariant under arbitrary shifts of f{z) by a real number, and thus we can choose 
its real part as we wish at a particularly convenient point. 

4 Properties of the critical solutions in the elliptic 
case and their Choptuik exponents 

We briefly present some of our results in the study of the critical solution in the elliptic 
case. The analysis follows closely the arguments in [3l H]. In all three non-linear systems 
(12-14) we have five singular points, z = ±0 represents the axis r = and regularity is 
imposed. The point z = oo represents the surface t = 0. Away from r = there is nothing 
special in this surface, hence regularity is also imposed. The simplest way to see that 
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^+ l/(o+)| |/(^+)| 



1.176 


2.609 


.892 


.364 


1.297 


2.674 


.903 


.392 


1.469 


2.771 


.910 


.397 


1.610 


2.822 


.914 


.404 


1.721 


2.871 


.918 


.413 


1.791 


2.929 


.923 


.427 


1.852 


2.998 


.928 


.439 



Table 1: Parameters determining the critical solution for dimensions 4 — 10, with dimension 
increasing from top to bottom. 

there is no problem at ^ = cxo is through a change of variables and a redefinition of the 
fields f{z),b{z) [3\. The singularities b{z±) = ±z± are those surfaces where the homothetic 
Killing vector becomes null. They correspond to the backward (resp. forward) light cone 
of the spacetime origin. For = z^ the solution should be smooth across this surface. 

The forward cone b{z-) = —Z- of the singularity represents the Cauchy horizon and we 
should not require more than continuity of /, b across this surface. We require smoothness 
of the spacetime below the forward cone and then extend the space by continuity inside it. 
In the elliptic case it is convenient to write f{z) = fm{z)e^^°-'^^\ requiring regularity both 
at the origin and at z^ leaves only four parameters to be determined: |/(0)|, w, z+, 
To find them we proceed by integrating from close to the origin towards z+ and also from 
a little below 2+ towards the origin, and match the results at an intermediate point, say 
z = 1. By matching the functions and the relevant first derivatives we can determine the 
unknown parameters in the critical solution. We do this for every dimension between d = 4 
and d = 10. The results are shown in Table 1. We did not optimise the numerical precision 
for these computations. In determining the analogous values of f{zJ) more precision is 
required because the solution is quite flat inside the forward light cone. The results will 
appear in fT6]. 

We compute the Choptuik exponents following the methods in [171 [HI [21 [!]• Given 
the critical solution, we can perturb it to find the critical exponent 7 as follows. Let h be 



6 



dimension 
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4 
5 
6 
7 
8 
9 
10 



261 
278 
287 
311 
352 
373 
.347 



Table 2: Choptuik exponents in dimensions 4—10 for the elliptic case. 



any function determining the critical solution: 6 or /. Next consider small perturbations 
around the critical solution: 



where hss{z) is the critical solution, e is a small number, /t is a constant, and /ipert(^) 
depends only on z. Substituting h{z, t) in the full equations of motion of motion (2,3) 
and keeping first order terms in e, gives a set of linear equations for the perturbation and 
an eigenvalue equation for n. This eigenvalue equation can in principle have a number of 
possible solutions for k. The solution with the largest value of Re{K) will be responsible 
for the fastest growing perturbation in (IT^ . and is called the "most relevant mode." The 
critical exponent is given by 



We require also that the perturbations are smooth at the singularities of the original equa- 
tions. This determines the hpc,j.t{z) up to some rescalings (which is fine given that the 
equations are linear) and also the most relevant value for Re(K). In fact, we take k, to be 
real. The values we obtain for the Coptuik exponents in dimensions from 4 — 10 appear 
in Table 2. The numerical accuracy decreases as the dimension increases. The third sig- 
nificant digit in Table 2 should not be trusted, specially in dimensions 9 and 10. More 
powerful numerical methods are necessary to get more accurate values for the exponents. 
We have preliminary results in the hyperbolic and parabolic cases indicating that there are 
also critical solutions in diverse dimensions with Choptuik exponents different from those 
found for the elliptic case [16]. 

As mentioned in the introduction the SL{2, R) symmetry of the classical type JIB string 
theory is supposed to break to SL{2, Z) once quantum effects are taken into account. This 
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raises the very interesting possibility that the critical solution in the quantum case will not 
be continuous self-similar but rather it will have discrete self-similarity as for the massless 
scalar field first analyzed by Choptuik jpLj. In other words, we want to know if there are 
elements T G SL{2,Z), such that: 

^ T ^ (I I) e SH2,Zi, (17) 

with Ar the corresponding echo parameter. To settle this question one would have to do the 
full numerical integration of Einstein's equation without assuming continuous self-similarity. 
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